clear;
m = 256;
n = 128;
% 256x128矩阵，每个元素服从Gauss随机分布
A = randn(m, n);
% 128x1矩阵，每个元素服从Gauss随机分布
% b = A * u, 则 Au - b = 0，u 即为我们要求的精确解
u = randn(n, 1);
b = A * u;

alpha = OLE_linprog(A, b);

% 绘制原信号
plot(u, 'b', 'LineWidth', 2);
hold on;
% 绘制恢复信号
plot(alpha, 'r.-');
grid on;
hold off;  
legend('Original', 'Recovery');
fprintf('\n恢复残差：');  
% 恢复残差 
norm(alpha - u, 2)